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Abstract. Consider the continuum of points along the edges of a network, 
i.e., a connected, undirected graph with positive edge weights. We measure 
the distance between these points in terms of the weighted shortest path 
distance, called the network distance. Within this metric space, we study 
farthest points and farthest distances. We introduce a data structure 
supporting queries for the farthest distance and the farthest points on two- 
terminal series-parallel networks. This data structure supports farthest- 
point queries in 0(k + logn) time after 0(n log p) construction time, 
where k is the number of farthest points, n is the size of the network, 
and p parallel operations are required to generate the network. 


1 Introduction 

Consider a geometric network with positive edge weights. For any two points 
on this network (i.e., points that may be vertices or in the interior of an edge), 
their network distance is the weight of a weighted shortest-path connecting them. 
Within this metric space, we study farthest points and farthest distances. We 
introduce a data structure supporting queries for the farthest distance and the 
farthest points on two-terminal series-parallel networks. 

As a prototype application, imagine the task to find the ideal location for a 
new hospital within the network formed by the streets of a city. One criterion for 
this optimization would be the emergency unit response time, i.e., the worst-case 
time an emergency crew needs to drive from the hospital to the site of an accident. 
However, a location might be optimal in terms of emergency unit response time, 
but unacceptable with respect to another criterion such as construction costs. We 
provide a data structure that would allow a decision maker to quickly compare 
various locations in terms of emergency unit response time. 

We obtain our data structure for two-terminal series-parallel networks by 
studying simpler networks reflecting parallel structure (parallel-path) and serial 
structure (bead-chains). Combining these insights, we support queries on flat 
series-parallel networks (abacus). Finally, we decompose series-parallel networks 
into a tree of nested abaci and combine their associated data structures. 

Our focus on supporting human decision makers with data structures deviates 
from the common one-shot optimization problems in location analysis, where 
we assume that only one factor determines suitable locations for some facility 
in a network. Moreover, we illustrate new ways of exploiting different parallel 
structures of networks that may be useful for tackling related problems. 



1.1 Preliminaries 


A network is defined as a simple, connected, undirected graph N = (V , E) with 
positive edge weights. We write w uv to denote the weight of the edge uv £ E 
that connects the vertices u,v £ V. A point p on edge uv subdivides uv into two 
sub-edges up and pv with w up = A w uv and w pv = (1 — A )w uv , for some A € [0, l]. 3 
We write p £ uv when p is on edge uv and p £ N when p is on some edge of N. 
The network distance between p,q £ N, denoted by d/v(p, q), is measured as 
the weighted length of a shortest path from p to q in N. We denote the farthest 
distance from p by diy(p ), i.e., diy{p) = max g£ jv d;sr(p,q )• Accordingly, we say a 
point p on N is farthest from p if and only if d]y{p,p) = djv(p). 

We develop data structures supporting the following queries in a network N. 
Given a point p on TV, what is the farthest distance from pi What are the farthest 
points from p in N1 We refer to the former as farthest-distance query and to the 
latter as farthest-point query. The query point p is represented by the edge uv 
containing p together with the value A £ [0,1] such that w up = A w uv . 



Fig. 1 . The operations (a) that generate two-terminal series-parallel networks (b). 

The term series-parallel network stems from the following two operations 
that are illustrated in Fig. 1. The series operation splits an existing edge uv into 
two new edges ux and xv where a: is a new vertex. The parallel operation creates 
a copy of an existing edge. A network N is two-terminal series-parallel when 
its underlying graph 4 can be generated from a single edge uv using a sequence 
of series and parallel operations; the vertices u and v are called terminals of N. 
We refer to the number of parallel operations required to generate N as the 
parallelism of N and to the number of series operations as the serialism of N. 

A network is called series-parallel when every bi-connected component is 
two-terminal series-parallel with respect to any two vertices. In this work, we 
only consider bi-connected networks; in future work, we shall adapt our treatment 
of multiple bi-connected components from cacti [5] to series-parallel networks. 

3 Observe that p V when A £ (0,1) in which case none of the sub-edges up and pv 
are edges in E. When A = 0 or A = 1, the point p coincides with u and v, respectively. 

4 The final graph is simple even if intermediate graphs have loops and multiple edges. 
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1.2 Related Work 


Duffin [8] studies series-parallel networks to compute the resistance of circuit 
boards. He characterizes three equivalent definitions of series-parallel networks 
and establishes their planarity. Two-terminal series-parallel networks admit linear 
time solutions for several problems that are NP-hard on general networks [3,15]. 
Since series-parallel networks have tree-with two [6], this applies to all problems 
with efficient algorithms on networks with bounded treewidth [1], 

A network Voronoi diagram subdivides a network depending on which site is 
closest [11] or farthest [9,14] among a finite set of sites. Any data structure for 
farthest-point queries on a network represents a network farthest-point Voronoi 
diagram where all points on the network are considered sites [4]. 

A continuous absolute center is a point on a network with minimum farthest- 
distance. Computing a continuous absolute center takes 0(n ) time on cacti [2] 
and 0(m 2 log n) time on general networks [12]. As a by-product, we obtain all 
continuous absolute centers of a series-parallel network in 0(n log p) time. 


1.3 Structure and Results 

We introduce a data structure supporting queries for the farthest distance and 
the farthest points on two-terminal series-parallel networks. We obtain this data 
structure by isolating sub-structures of series-parallel networks: In Sections 2 
and 3, we study networks consisting of parallel paths and networks consisting of 
a cycle with attached paths (bead-chains), respectively. In Section 4, we combine 
these results into abacus networks, which are series-parallel networks without 
nested structures. Finally, we combine these intermediate data structures to 
obtain our main result in Section 5. Table 1 summarizes the characteristics of 
the proposed data structures and compares them to previous results. 


Type Farthest-Point Query 

Size 

Construction Time Reference 

General 

0(k + logn) 

0(m 2 ) 

0(m 2 logn) 

[4] 

Tree 

0(k) 

O(n) 

O(n) 

[5] 

Cycle 

O(logn) 

O(n) 

0(n) 

[5] 

Uni-Cyclic 

0(k + logn) 

O(n) 

0(n) 

[5] 

Cactus 

0(k + logn) 

O(n) 

O(n) 

[5] 

Parallel-Path 

0(k + logn) 

O(n) 

0(n ) 

this work 

Bead-Chain 

0(k + logn) 

O(n) 

O(n) 

this work 

Abacus 

0(k + logn) 

O(n) 

0(n log p) 

this work 

Series-Parallel 

0(k + logn) 

O(n) 

0(n log p) 

this work 


Table 1. The traits of our data structures for queries in different types of networks, 
with n vertices, m edges, k reported farthest points, and parallelism p. 
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2 Parallel-Path Networks 


A parallel-path network consists of a 
bundle of edge disjoint paths connect¬ 
ing two vertices u and v, as illustrated 
in Fig. 2. In terms of series-parallel 
networks, parallel-path networks are 
generated from an edge uv using a se¬ 
quence of parallel operations followed 
by a sequence of series operations. 

Let Pi,P 2 ,..., P p be the paths of 
weighted lengths w\ < W 2 < • • • < w p 
between the terminals u and v in a parallel-path network N. Consider a shortest 
path tree 5 from a query point q G N. As depicted in Fig. 3, there are three 
cases: either all shortest paths from q reach v via u ( left case), or all shortest 
paths from q reach u via v ( right case), or neither ( middle case). We distinguish 
the three cases using the following notation. Let Xi denote the farthest point 
from x € {it, u} among the points of path Pi, i.e., Xi is a point on Pi such that 
d(x,Xi) = max ye p i d(x,y). Together with Fig. 3, the next lemma justifies our 
choice of the names left case, middle case, and right case. 


ul _ Jr 

V jTdi 

<Ti S3. rfC 

- - -'Ti, 

(a) The left case. 

Fig. 3. The three cases for queries in parallel-path networks. Consider the shortest path 
tree from a query point q € Pi along the paths Pi (center), Pi (bottom), and P, (top). 
In the left case (a), we reach v via u. In the right case (c), we reach u via v. In the 
middle case (b), neither holds, i.e., we enter path Pi from both terminals u and v. 
Points colored red are reached fastest via a path through u or towards u along uv, while 
points colored blue are reached fastest via a path through v or towards v along uv. 

Lemma 1. Consider a query q from the i-th path of a parallel-path network. 

(i) We are in the left case when q lies on the sub-path from u to v j with q ^ Vi, 

(ii) we are in the middle case when q lies on the sub-path from to Ui, and 
(Hi) we are in the right case when q lies on the sub-path from Ui to v with q ^ Ui. 

5 More precisely, we consider extended shortest path trees [14] which result from splitting 
each non-tree edge st of a shortest path tree into two sub-edges sx and xt, where all 
points on sx reach the root through s and all points on xt reach the root through t. 




Fig. 2. A parallel-path two-terminal series- 
parallel network with parallelism p = 7. 
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Proof. Assume we are in the left case, where every shortest path from the query 
point q £ Pi to v contains u, i.e., d(q,v) = d{q,u) + d{u,v) and d{q,v) < w qv . 
The latter implies q ^ Vj, since d(v,v j) = Wy iV . Moreover, v, cannot lie on the 
sub-path from q to u along Pi, since otherwise 

_ _ q^Vi _ 

d(q,v) = d(q,Vi) + d(vi,v) > d(vi,v ) , 

contradicting the choice of V{ as farthest point from v on P t . Therefore, q must 
lie between u and v, along Pi when we are in the left case. 

Conversely, assume q lies between u and Vi along Pi with q ^ tij. No shortest 
path from q to v can contain v t in its interior. Hence, every shortest path from 
q to v reaches v via u, i.e., the left case applies. Symmetrically, the right case 
applies if and only if q ^ Ui lies between Hi and v. Consequently, only the middle 
case remains for all points on the sub-path from v l to u t and the claim follows. □ 

Using Lemma 1, we deal with the three cases as follows. 


Left Case and Right Case In the left case, every shortest path from q £ Pi to 
any point outside of Pi leaves P,; through u. Hence, the farthest point from q on 
Pj with j ^ i is the farthest point Uj from u on Pj. The distance from q to qj is 
d(q , qj) = d(q, u) + d(u, Uj) = w qu + On the other hand, the farthest point 

qi from q on Pi itself moves from u, to v as q moves from u to Vi maintaining a 
distance of d{q,qi) = Wl + Wi . Therefore, the farthest distance from q in N is 


d(q) = max 


W i + Wi 


2 

Wi + w p 
2 


, max w qu + 
j 


Wi + Wj 


iii^p 


= < 


v qu 


Wi+Wp-I Wp-Wp -1 

-I-—— — if i = p and ————— 


„ _ » ana ———-— < w au 

2 2 ~ q 


w i + w p 


if i = p and 


w p - w p -1 


> w qu 


The first case means that, for queries from anywhere other than P p , the farthest 
points lie on the u-t>-paths of maximum length. The second and third case distin¬ 
guish whether P p contains a farthest point for queries from P p itself. Accordingly, 
we answer a farthest point query from q £ Pi in the left case as follows. 

— If i ^ p, we report all Uj where Wj = w p and i ^ j. 

— If i = p and Wp f’’- 1 < Wq U , we report all Uj where Wj = w p -\ and j ^ p. 

— If i = p and Wp f' p ^ 1 > w qu , we report the farthest point q p from q on P p 
using a binary search along the sub-path of P p from u p to v. 

The overlap between the last two cases covers the boundary case when a query 
from P p yields a farthest point on P p itself and farthest points on other u-u-paths. 

Swapping u and v above yields the procedure for the right case. Thus, an¬ 
swering farthest-point queries takes 0(k + logn) time in the left and right cases. 
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Middle Case In the middle case, there are no farthest points from q on Pi itself 
and every path Pj with j ^ i contains points that we reach from q via u as well 
as points that we reach from q via v. Let qj be the farthest point from q along the 
cycle formed by Pj and Pj. Since the distance from q to qj is d(q,qj) = W, + Wj ; 
the farthest distance d(q) from q in TV is 


d(q) = max 


! Wi + w p 
~ 

Wp + Wp -1 
2 


if i ^ p 
if i = p 


The first case applies for queries from anywhere other than P p who have their 
farthest points on the longest u-u-paths, i.e., on the paths Pj with Wj = w p . 
The second case applies for queries from P p who have their farthest points on 
the second longest u-u-paths, i.e., on the paths Pj with Wj = w p -\. Using binary 
search, we can answer a farthest point query from q £ Pi in the middle case by 
reporting the points qj on those k paths Pj that contain farthest points from q. 
To improve the resulting query time of 0(k log n) , we take a closer look at the 
position of qj relative to Uj and Vj. As illustrated in Fig. 4, the farthest point qj 
from q £ Pi along Pj moves from Uj to Vj as q moves from Vj to Uj. 


Lemma 2. Let N be a parallel-path network with terminals u and v. For any 
point x £ N, let Xi denote the farthest point from x along the u-v-path Pj. 

(i) The sub-path from Vi to Ui has length d(u,v). 

(ii) For every point q along the sub-path from Vj to Uj, the sub-path from Vj to q 
has the same length as the sub-path from Uj to qj for any j ^ i. 

Proof. We have w U y i = w V u t for every path Pj, since 

d(u, v ) + Wj 

WuVi T Wy i u i - W U u i = ~ = W V y i — Wyy i T Wji i y i . 

The above implies the first claim, since adding a nutritious zero yields 


Wy 


= Wy 


i+iWuVi-Wytii) = d(u,Ui)-W vni = d(u,v)+w V u i -w V u i =d(u,v) 


For the second claim, we compare the length of the shortest path from q to 
qj via v with the one via u. By studying Fig. 4, we obtain the following. 

d{q, Qj) = Wqy, + Wy iU + W U y. + Wy jUj ~ Wqj Uj (1) 

d{q, qj) = Wqjijj + w V uj + w vUi + Wy iUi - w q Ci (2) 

We generate nutritious zeros from the identities w U y i = w v a i , w U y . = w v s., and 
Wu i y i = d(u,v) = Wjj.y .. Rearranging terms yields the second claim, as 

=o =o =o =o 

‘Z'Wqvi = ‘^'Wqvi ”1“ ^VjU ^vui ^uijj Wvuj H - WVjUj WVjUi H ”qjUj ^qjUj ) 

— d(q, qj) - [Wq jUj + w vUj + w vUi + Wy iUi ~ w q y z ] + 2 Wq jUj 
( 2 ) 

= d(q, qj) - d(q , q 3 ) + 2 Wq j u j = 2 w qjUj • □ 
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Using Lemma 2, we interpret the searches 
for ft on the sub-path from ft to Uj, as a 
single search with a common key q in mul¬ 
tiple lists (the ft-ft-sub-paths) of compara¬ 
ble search keys (the vertices along these sub¬ 
paths). Using 0(n) time, we construct a frac¬ 
tional cascading data structure [7] supporting 
predecessor queries on the sub-paths from ft 
to Uj for those paths Pj where w 3 = w p - 1 . 6 

We answer a farthest-point query from 
q € Pi as follows. If i ft p, we locate and 
report q p along P p in O(logn) time. If i = p 
or w p = w p - 1 , the remaining farthest points from q are the qj where j / j and 
Wj = w p - 1 ; we report them in 0{k logn) time using the fractional cascading 
data structure. This query might report a point on Pi, which would be ft for 
queries from outside P,;. For queries from within Pi, we omit this artifact. 



Fig. 4. The positions of the points 
along the cycle Pi 


Theorem 1. Let N be a parallel-path network with n vertices. There is a data 
structure with 0(n) size and 0(n ) construction time supporting 0(k + logn )-time 
farthest-point queries on N, where k is the number of farthest points. 


3 Bead-Chain Networks 

A bead-chain network consists of a main cycle 
with attached arcs so that each arc returns 
to the cycle before the next one begins. An 
example is depicted in Fig. 5. Bead-chains are 
series-parallel networks where we first subdi¬ 
vide a cycle using series operations, then we 
apply at most one parallel operation to each 
edge of this cycle followed by series operations 
that further subdivide the arcs and cycle. 

Consider a bead-chain network N with 
main cycle C and arcs oq,..., a s . Let aq and 
bi be the vertices connecting C with the i-th 
arc. Without loss of generality, the path ft 
from ai to bi along C is at most as long as at. 

Otherwise, we swap the roles of ai and ft. 

We first study the shape of the function ft (a;) that describes the farthest 
distance from points along the main cycle to any point on the Pth arc, i.e., 
ft (a;) = maxygc,. d{x,y). When considering only the i-th arc, we have a parallel- 
path network with three paths. Let x denote the farthest point from x £ C on C 
itself and let ft denote the farthest point from x on arc ai. From the analysis in 
the previous section, we know that ft (a;) has the shape depicted in Fig. 6. 

6 We consider paths of length w v -\ instead of w p , because we treat P p separately. 



Fig. 5. A bead-chain with four arcs 
(colored) around its cycle (black). 
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d(ai,bi) 



Fig. 6. The shape of the function di(x) describing the distance from x € C to the 
farthest point Xi from x among the points on the i-th arc at. 


When walking along the main cycle, we encounter ai, bi, di, and bi in this 
order or its reverse. From a,; to bi, the point x t moves from dj to bi maintaining 
a constant distance. From bi to di, the point £7 stays at bi increasing in distance. 
From di to bi, the point Xi moves from bi back to di, again, at a constant distance. 
Finally, 27 stays at a, with decreasing distance when x moves from bi to ai. 

Since the farthest distance changes at the same rate when we move towards 
or away from the current farthest-point, the increasing and decreasing segments 
of any two functions di and dj have the same slope except for their sign. 

The height of the upper envelope D of the functions di,... ,d s at q € C 
indicates the farthest distance from q to any point on the arcs and the i-th arc 
contains a farthest point from q when di coincides with D at q. We construct D 
in linear time using the shape of the functions d\,..., d s described above. 

We need to consider an arc separately from the other arcs when it is too long. 
We call an arc on overlong when the path /3j is longer then remainder 7 \ of the 
main cycle. Figure 7 illustrates an overlong arc on with its function di. 



d(a,i,bi) 
d(cii, Cti ) 


a i 5 . di bi a i 

(b) The farthest distance to an overlong arc. 



Fig. 7. An overlong arc ai (blue) in a bead-chain network where /3i (green) is longer 
then the remaining cycle 7 * (orange). The shape of di is the same as for non-overlong 
arcs, but its high plateau may horizontally overlap with the high plateau of other arcs. 












Lemma 3. Every bead-chain network N has at most one overlong arc and the 
functions d\,..., d s of the remaining arcs oq ,a s satisfy the following. 

(i) The high plateaus of d \,..., d s appear in the order as their arcs oq,..., a s 
appear along the cycle and no two high plateaus overlap horizontally. 

(ii) The low plateaus of d\,... ,d s appear in the order as their arcs oq ,... ,a s 
appear along the cycle and no two low plateaus overlap horizontally. 

Proof. The existence of two overlong arcs aq and ak leads to the contradiction 
wp i > w li > wp k > w lk > wp i , since 7 j contains /3k and 7 *, contains fa. 

Let ot \,..., a s be the non-overlong arcs of N as they appear along the cycle. 
For every arc oq, with i = 1,..., s, the farthest points ai and bi from the endpoints 
a,; and bi of oq appear along 7 i. Therefore, the points di, 61 , a^, fa, ■ ■ ■, d s , and b s 
appear in this order along the cycle C. Claim (i) follows, since the high plateau 
of di lies between di and Claim (ii) follows, since any non-overlong arc oq has 
its low plateau along fa and since fa ,..., fa are separate by definition. □ 

As suggested by Lemma 3, we incrementally construct the upper envelope of 
the functions di corresponding to non-overlong arcs and treat a potential overlong 
arc separately. When performing a farthest-point query from the cycle, we first 
determine the farthest distance to the overlong arc and the farthest distance to 
all other arcs. Depending on the answer we report farthest points accordingly. 

Lemma 4. Let aq,..., a s be the arcs of a bead-chain that has no overlong arc. 
Computing the upper envelope D of d\,... ,d s takes O(s) time. 

Proof. We proceed in two passes: in the first pass, we consider only the high 
plateaus and non-constant segments of d \,..., d s ; the respective low plateaus are 
replaced by extending the corresponding non-constant segments. In the second 
pass, we traverse the partial upper envelope from the first pass again and compare 
it with the previously omitted low plateaus, thereby constructing D. 

Let d/ be the function resulting from replacing the low plateau of di by 
extending its non-constant segments and let D' i be the upper envelope of d[,..., df 
In the first pass, we construct D' s incrementally. Assume we have D' i _ 1 and would 
like to obtain D[ by inserting d\ into Dfa 1 , as illustrated in Fig. 8 . We perform 
this insertion by walking from di —the left endpoint of the high plateau of dj- in 
both directions updating the current upper envelope D t _^. Locating a,; takes 
constant time, since di is the first bending point to the right of 6 ,_ 1 . 

There is no more than one increasing segment of D' i _ l between di and bi. 
Assume, for a contradiction that there are two increasing segments Si and S 2 
between a,; and bi. Neither of them has their higher endpoint between a,; and bi, 
since there would be two horizontally overlapping high plateaus, otherwise. Since 
si and S 2 have the same slope, only one of them can be part of the upper envelope.' 
Analogously, the same holds for decreasing segments. Therefore, inserting the 
high plateau of di into the upper envelope D' t i takes constant time. 

' When the segments si and S 2 happen to overlap, we consider only the segment of 
the arc distance function that was inserted first to be part of the upper envelope. 
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Fig. 8. An incremental step where we construct D[ from D [_i (black) and d[ (orange). 
The treatment of the low plateau of di (dashed) is deferred to the second pass. 


If the decreasing segment of d[ appears along D[ at all, then it appears at bi. 
We update the previous upper envelope D' i _ 1 by walking from bi towards oq until 
the decreasing segment of d' vanishes below D[_ 1 . We charge the costs for this 
walk to the segments that are removed from the previous upper envelope D' i _ 1 . 
We proceed in the same fashion with the increasing segment of d! r by walking 
from di towards bi. Each non-constant segment appears at most once along any 
intermediate upper envelope and is never considered again after its removal. 
Therefore, the total cost for inserting all non-constant segments—and, thus, the 
total cost of constructing D ' s —amounts to O(s). 

In the second pass, we construct the desired upper envelope D from D' s . Since 
no two low plateaus overlap horizontally, we simply walk along D' s comparing 
its height to the height of the current low plateau, if any. This takes O(s) time, 
since D' s has O(s) bending points and since there are s low plateaus. □ 

To answer a farthest-point query from q £ C, we need to find its farthest arcs, 
i.e., the arcs containing farthest points from q. Suppose each point along the main 
cycle C has exactly one farthest arc. Then we could subdivide C depending on 
which arc is farthest and answer a farthest-arc query by identifying the function 
among d\,.... d s that defines the upper envelope D on the sub-edge containing q. 
On the other hand, there could be multiple farthest arcs when several functions 
among d±,... ,d s have overlapping increasing or decreasing segments. In this case, 
we could store the at most two farthest arcs from plateaus directly with the 
corresponding segments of D. However, storing the farthest arcs from increasing 
and decreasing segments directly would lead to a quadratic construction time. 
Instead, we rely on the following observation. An arc a is considered relevant 
when there exists some point x £ C such that a is a farthest arc for x and a is 
considered irrelevant when there is no such point on the main cycle. 
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Lemma 5. Let ai, aj, and ak be arcs that appear in this order in a bead-chain 
without overlong arc. The arc aj is irrelevant when ct^ and ak are farthest arcs 
from some query point q such that di and dk are both decreasing/increasing at q. 

Proof. As illustrated in Fig. 9, q lies between at and bk, since di and dk are both 
decreasing at q. Hence, q lies between aj and bj , i.e., dj is decreasing at q, as well. 
We show dj(x) < dfx) for all x £ C, which implies that aj is irrelevant. 



Fig. 9. The constellation from Lemma 5, where ai and a*, are farthest arcs from q where 
the farthest distances decreases as q moves in clockwise direction (a). A comparison of 
the arc distance functions di and dj reveals that aj is irrelevant in this case (b). 


Consider the difference A{x) := di{x) — dj(x). We have A(a,i) = A(q) > 0, 
as di and dj are decreasing with the same slope from q to aj. We observe that 
A{bi) = A(ai) + d(ai,bi), as di remains constant from ai to 6 j while dj decreases. 
We have A(aj) = A(bi)+2d(bi, aj), as di increases from to aj while dj decreases. 
We continue in this fashion and obtain the following description of A. 


A(ai ) = A(q) 

A(bi) = A(ai) + d(ai , h) 
A(aj) = A(bi ) + 2 d(bi,aj) 
A(bj) = A(aj) + d{aj,bj) 


A{di) = A(bj) 

A(bi) = A(di) - d(di,bi) 
A(dj) = A(bi) - 2 d(bi, dj) 
A(bj) = A(bj) — d(dj,bj) 


The claim follows as the above implies A(x) > A(q) > 0 for all x £ C, due to the 
symmetries d(ai , bf) = d(di, bf), d(bi, aj) = d(bi , dj), and d(aj,bj) = d(dj,bj). □ 


Corollary 1. Let q be a point on the cycle of a bead-chain with no overlong arc. 
The farthest arcs from q that correspond to decreasing/increasing segments of D 
form one consecutive sub-list of the circular list of relevant arcs. 
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Using Corollary 1, we answer farthest-arc queries from the main cycle of a 
bead-chain network without overlong arc as follows. When D has a plateau at the 
query point g, we report the at most two farthest arcs stored with this plateau. 
When D has an increasing/decreasing segment at g, we first report the farthest 
arc a that is stored directly with this segment. We report the remaining farthest 
arcs by cycling through the circular list of relevant arcs starting from a in both 
directions until we reach a relevant arc that is no longer farthest from q. 

Theorem 2. Let TV be a bead-chain network with n vertices. There is a data 
structure with 0(n ) size and 0(n) construction time supporting 0(k + log n)-time 
farthest-point queries on TV, where k is the number of farthest points. 

Proof. Let TV be a bead-chain network with main cycle C and arcs a, aq,..., a s 
where only a may be overlong. Let N\a be the network obtained by the removing 
the potentially overlong arc from TV, i.e., TV \ a := C U a\ U • • • U a s . 

We support queries from the main cycle C by constructing (i) a data structure 
for queries in the parallel-path network a U C and (ii) a data structure for queries 
in the bead-chain network N \ a that has no overlong arc. This construction 
takes linear time, due to Theorem 1 and Lemma 4. For a farthest-point query 
from q £ C, we perform a farthest-distance query from q in a U C and in TV \ a. 
Depending on which of these two queries reported the largest farthest distance, 
we conduct the appropriate farthest-point queries in a U C or in TV \ a. 

We support queries from arc cq (from a) by combining (i) a data structure for 
farthest-point queries in the cycle cq U /3j (in a U C), with (ii) the data structure 
supporting farthest-point queries in TV from the cycle C described above, and 
(iii) the partial upper envelope D' from the proof of Lemma 4. 8 

We begin a query from q £ cti with the corresponding query from q in a,; U /?,. 
After that, we report any farthest points from q in TV outside of a,; based on the 
position of q along cq: When q lies between a.j and bi (left case), we report the 
farthest points from q in TV \ cq with a query from a* in TV. When q lies between 
hi and bi (right case), we report the farthest points from q in TV \ a,; with a 
query from bi in TV. When q lies between bi and hi (middle case), we attempt a 
query from the unique point q' on /?.; with d(bi,q) = d(ai,q). Since q' preserves 
the relative position from q on aq, the query from q' yields the farthest points 
from q in TV outside of a,—except when q' has only a single farthest point on aq. 
This occurs when D has the low plateau of di at q'. Fortunately, we are able to 
recover the correct answer, since the upper envelope of d ±,..., dj_i, d*+i,..., d s 
coincides at q' with the partial upper envelope D' from the proof of Lemma 4. 

Our data structure has 0(n) size and construction time, since each vertex 
appears only in a constant number of sub-structures each of which have linear 
construction time. Every farthest-point query takes 0(/c + logn) time, because 
each query consists of a constant number of 0(log?r)-time farthest-distance 
queries followed by a constant number of farthest point queries only in those 
sub-structures that actually contain farthest points from the original query. □ 

8 Recall that D' was the upper envelope of the functions that result from replacing 
the low plateaus of di, ..., d s by extending their increasing/decreasing segments. 


12 



4 Abacus Networks 


An abacus is a network A consisting of a parallel-path network N with arcs 
attached to its parallel paths, as illustrated in Fig. 10. Let P\, ... , P p be the 
parallel paths of N and let Bi be the i-th parallel path with attached arcs. 



Fig. 10. An abacus with the arcs (colored) attached to its parallel-path network (black). 

We split farthest-point queries in an abacus into an inward query and an 
outward query, an inward query considers farthest points on the bead-chain 
containing the query point; an outward query considers farthest points on the 
remaining bead-chains. We first perform the farthest distance version of inward 
and outward queries before reporting farthest points where appropriate. Figure 11 
illustrates how we treat inward and outward queries in the following. 

For an inward query from q on Bi , we construct the bead-chain network B[ 
consisting of Bi with an additional edge from u to v of weight d(u, v), as illustrated 
in Fig. 11a. Since B[ preserves distances from A, the farthest points from q € B' 
are the farthest points from q among the points on Bi in A. 

For outward queries in an abacus, we distinguish the same three cases as 
for parallel-path networks: we are in the left case when every shortest path tree 
reaches u before v, we are in the right case when every shortest path tree reaches 
v before u, and we are in the middle case otherwise. Analogously to Lemma 1, 
the left case applies when we are within distance d(u,Vi) from u and the right 
case applies when we are within distance d(v,Ui) from v. 

For an outward query from q € Bi in the left case, q has the same farthest 
points as u outside of Bi. During the construction of the networks B [,..., B' p 
for inward queries, we determine a list L ? of the farthest points from u in Bj. 
Similarly to our treatment of the left case for parallel-path networks, we only 
keep the list achieving the highest farthest distance and the lists achieving the 
second highest farthest distance. With this preparation, answering the query 
for q amounts to reporting the entries of the appropriate lists Lj with j i. 
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(b) The left case of an outward query. 




(c) The right case of an outward query. (d) An outward query from an arc. 




(f) Collapsing the chains. 


Fig. 11. Inward (a) and outward (b-f) queries for the abacus network from Fig. 10. 
Inward queries are answered in the bead-chain containing the query (a). Outward 
queries in the side case are answered with queries form the terminals (b,c). Outward 
queries in the middle case from arcs are translated to queries from the path (d) and 
then to queries from a virtual edge (e). From the perspective of the virtual edge, we 
conceptually collapse all bead-chains of the abacus to support virtual queries (f). 
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For middle case outward queries, we proceed along the following four steps: 
First, we translate every outward query from an arc of Bj to an outward query 
from the path P t . i.e., we argue that it suffices to consider outward queries from 
the parallel paths (Fig. lid). Second, we translate outward queries from T\ to 
outward queries from a virtual edge e connecting the terminals (Fig. lie). Third, 
we speed up queries from the virtual edge by superimposing the data structures 
for the bead-chains B\ Ue, ..., B p Ue, i.e., by conceptually collapsing the parallel 
chains (Fig. Ilf). Finally, we recover the correct answer to the original outward 
ciuery from the answer obtained with an outward query from the virtual edge. 

Lemma 6. Let a be an arc in an abacus and let (3 be the other path connecting 
the endpoints of a. For every point q € a in the middle case, there is a point 
q' £ (3 such that q' has the same outward farthest points as q. 

Proof. We continuously deform a to f3 maintaining the relative position of q to 
the endpoints of a. The distance from q to all points in the network decreases at 
the same rate, hence, the outward farthest points remain the same. □ 

We introduce a virtual edge e from u to v of length w p , i.e., the length of 
the longest u-u-path P p in the underlying parallel-path network, as illustrated 
in Fig. lie. Let u be the farthest point from u on e and let v be the farthest 
point from v on e. From Lemma 2, we know that the sub-edge uv of e has length 
d(u, v) and, thus, the same length as the sub-path from Hi to Vi on each parallel 
path Pi. We translate an outward query from q £ Pi to a query from the unique 
point q on e such that q has the same distance to u and to v as q to Ui and to u,. 

Lemma 7. For q £ Pi in the middle case, the farthest points from q in Pi U Bj 
are the farthest points from q in e LI Bj for every j ^ i. 

Proof. We continuously elongate Pi to e maintaining the relative position of q to 
u and v and, thus, to iq and Hi. At the end of this process q coincides with q. 
The distance from q to all points outside of Bi increases uniformly. Hence, the 
outward farthest points remain the same throughout the deformation. □ 

It would be too inefficient to inspect each bead-chain network Bj U e with 
j i to answer an outward query from q £ Pi. Instead, we first determine 
the upper envelopes of the farthest-arc distances D i,..., D p along e in each 
B\ U e,..., B p U e and then compute their upper envelope U\ as well as their 
second level C/ 2 , i.e., the upper envelope of what remains when we remove the 
segments of the upper envelope. Computing the upper envelope and the second 
level takes O(nlogp) time, e.g., using plane sweep. Using fractional cascading, 
we support constant time jumps between corresponding segments of U\ and C/ 2 . 
The resulting structure occupies 0(n) space, since each of the 0(n) arcs along 
any bead-chain contributes at most four bending points to C /1 and C/ 2 . 

We answer an outward query from q £ Pi in the middle case by translating q 
to q. When the segment defining U\ at q is from some arc a of Bj with j ^ i, 
then a contains an outward farthest point from q. When the segment defining 
U\ at q corresponds to an arc of Bi, then we jump down to C/ 2 , which leads us 
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to an arc containing an outward farthest point from q. We report the remaining 
arcs with outward farthest points by walking q along Ui and U^- In order to skip 
long sequences of segments from we introduce pointers along Ui to the next 
segment from another chain in either direction. Answering outward queries in 
the middle case takes 0(k + logn) time after 0(nlogp) construction time. 

Theorem 3. Let N be an abacus with n vertices and p chains. There is a data 
structure of size 0(n) with 0(nlogp) construction time supporting farthest-point 
queries on N in 0(k + logn) time, where k is the number of farthest points. 


5 Two-Terminal Series-Parallel Networks 

Consider a two-terminal series-parallel network N. By undoing all possible series 
operations and all possible parallel operations in alternating rounds, we reduce N 
to an edge connecting its terminals and decompose N into paths that reflect its 
creation history. The colors in Fig. 12 illustrate this decomposition. 



Fig. 12. A two-terminal series-parallel network with colors indicating the parallel 
operations in a potential creation history: starting with a single red edge, we create a 
parallel yellow and a parallel blue edge. Then we subdivide the blue edge using series 
operations until we create a parallel purple edge for one of the blue edges and so forth. 


Lemma 8. Let N be a series-parallel network with parallelism p and serialism s. 
Identifying the terminals of N and reconstructing its creation takes 0(s+p) time. 

Proof. Series-parallel networks are planar [ 8 ]. We maintain a series-parallel net¬ 
work N together with its dual N* throughout the following reduction process 
illustrated in Fig. 13. We keep two arrays d and d* with n entries each where d[i\ 
stores a list of the vertices of N with degree i and d* [i] stores a list of the vertices 
of N* (faces of N) with degree i. Each vertex of N and each face of N maintain 
a pointers to their position in the lists to facilitate constant time deletions. 

We proceed in alternating rounds where we either reverse as many series 
operations or as many parallel operations as possible in each round. To reverse 
series operations we delete all degree two vertices of N. 


16 








(a) Step 1: reverting series operations. (b) Step 2: reverting parallel operations. 





(e) Step 5: reverting series operations. 


(f) Step 6: reverting parallel operations. 


Fig. 13. The reduction process for the two-terminal series-parallel network from Fig. 12. 
In each step, we contract all degree two vertices in the primal network to revert series 
operations or we contract all degree two vertices in the dual (blue) to revert parallel 
operations. When reverting parallel operations, the colors indicate how we backtrack 
the creation history in Fig. 12. For the sake of clarity, the dual vertex for the outer face 
(blue square) is only shown in the last step where no further reduction is possible. 
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Each deletion changes the degree of two faces of N that were incident to the 
removed vertex so we move these faces to their new positions in d*. To reverse 
parallel operations we proceed in the exact same fashion by removing all degree 
two vertices of N* and updating d accordingly. We recover the creation history of 
N by keeping track of when parallel edges were removed during the reversal of a 
parallel operation. This reduces N to the edge connecting its terminals in 0(s+p) 
steps, since reversing each of the s + p operations takes constant time. □ 

Once we know the terminals u and v of N, we compute the shortest path 
distances from u and from v in 0(n log p) time. 9 Consulting the creation history, 
we determine a maximal parallel-path sub-network P of N with terminals u and 
v. As illustrated in Figs. 14 and 15, every bi-connected component X of N that 
is attached to some path of P between vertices a and b is again a two-ternrinal 
series-parallel network with terminals a and b. We recurse on these bi-connected 
components. When this recursion returns, we know a longest a-&-path in X and 
attach an arc from a to b of this length to P. The resulting network is an abacus 
A. The abaci created during the recursion form a tree T with root A. Alongside 
with this decomposition we also create our data structures for the nested abaci. 

The size of the resulting data structure remains O(n), since the data structure 
for each nested abaci consumes space linear in the number of its vertices and 
each vertex in any nested abaci can be charged to one of the series or parallel 
operations required to generate the original network. 

We translate a query q to a query in the abacus A; queries from a bi-connected 
component X attached to P in TV will be placed on the corresponding arc of 
A. Whenever the query in A returns a farthest point on some arc, we cascade 
the query into the corresponding nested data structure. We add shortcuts to 
the abacus tree T in order to avoid cascading through too many levels of T 
without encountering farthest-points from the original query. This way, answering 
farthest-point queries in N takes 0(k + logn) time in total. 

Theorem 4. Let N be a two-terminal series-parallel network with n vertices and 
parallelism p. There is a data structure of size 0(n ) with 0(n\ogp) construction 
time that supports 0(k + logn)-time farthest-point queries from any point on N, 
where k is the number of farthest points. 

6 Conclusion and Future Work 

In previous work, we learned how to support fartlrest-point queries by exploiting 
the treelike structure of cactus networks [5]. In this work, we extended the arsenal 
by techniques for dealing with parallel structures, as well. In future work, we aim 
to tackle more types of networks such as planar networks, /c-almost trees [10], 
or generalized series-parallel networks [13]. Moreover, we are also interested in 
lower bounds on the construction time of data structures supporting efficient 
farthest-point queries to guide our search for optimal data structures. 

9 The priority queue in Dijkstra’s algorithm manages never more than p entries. 
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(a) A two-terminal series-parallel network. 



Fig. 14. Decomposing a two-terminal series-parallel network (a) into nested abaci. Each 
involved network is two-terminal series-parallel and we mark the terminals with empty 
circles. Replacing the nested structures in each involved network (b-g) yields the abacus 
network shown to its right. Each colored cycle consists of the shortest path and the 
longest path connecting the terminals of each nested structure X ; their weighted lengths 
determine the weight of the corresponding arc replacing X. This process terminates 
when we reach an abacus (d, g, i), i.e., when there are no more nested structures. 
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Fig. 15. The tree of nested abaci for the two-terminal series-parallel network from 
Fig. 14. The inner nodes of this tree correspond to two-terminal networks with nested 
structures that are indicated with colors; the leaves correspond to abacus networks 
without nested structures. A query would start at the root abacus and cascade into 
nested structures when necessary. For instance, when a query at the root abacus yields 
a farthest point on the blue arc and a farthest point on the red arc, we would perform 
subsequent queries in the abaci stored in the left and middle child of the root. 


20 









































































References 


1. Arnborg, S., Lagergren, J., Seese, D.: Easy problems for tree-decomposable graphs. 
Journal of Algorithms 12(2), 308-340 (1991) 

2. Ben-Moshe, B., Bhattacharya, B., Shi, Q., Tamir, A.: Efficient algorithms for center 
problems in cactus networks. Theoretical Computer Science 378(3), 237 - 252 (2007) 

3. Bern, M.W., Lawler, E.L., Wong, A.L.: Linear-time computation of optimal sub¬ 
graphs of decomposable graphs. Journal of Algorithms 8(2), 216-235 (1987) 

4. Bose, P., Dannies, K., De Carufel, J.L., Doell, C., Grimm, C., Maheshwari, A., 
Schirra, S., Smid, M.: Network farthest-point diagrams. Journal of Computational 
Geometry 4(1), 182-211 (2013) 

5. Bose, P., De Carufel, J.L., Grimm, C., Maheshwari, A., Smid, M.: Optimal data 
structures for farthest-point queries in cactus networks. Journal of Graph Algorithms 
and Applications 19(1), 11-41 (2015) 

6 . Brandstadt, A., Le, V.B., Spinrad, J.P.: Graph Classes: A Survey. SIAM Monographs 
on Discrete Mathematics and Applications, Society for Industrial and Applied 
Mathematics (SIAM), Philadelphia, PA (1999) 

7. Chazelle, B., Guibas, L.J.: Fractional cascading: I. A data structuring technique. 
Algorithmica 1(2), 133-162 (1986) 

8 . Duffin, R.J.: Topology of series-parallel networks. Journal of Mathematical Analysis 
and Applications 10(2), 303-318 (1965) 

9. Erwig, M.: The graph Voronoi diagram with applications. Networks 36(3), 156-163 

( 2000 ) 

10. Gurevich, Y., Stockmeyer, L.J., Vishkin, U.: Solving NP-hard problems on graphs 
that are almost trees and an application to facility location problems. Journal of 
the ACM 31(3), 459-473 (1984) 

11. Hakimi, S.L., Labbe, M., Schmeichel, E.: The Voronoi partition of a network and 
its implications in location theory. ORSA Journal on Computing 4(4), 412-417 
(1992) 

12. Hansen, P., Labbe, M., Nicolas, B.: The continuous center set of a network. Discrete 
Applied Mathematics 30(2-3), 181 195 (1991) 

13. Korneyenko, N.M.: Combinatorial algorithms on a class of graphs. Discrete Applied 
Mathematics 54(2-3), 215-217 (1994) 

14. Okabe, A., Satoh, T., Furuta, T., Suzuki, A., Okano, K.: Generalized network 
Voronoi diagrams: Concepts, computational methods, and applications. Interna¬ 
tional Journal of Geographical Information Science 22(9), 965-994 (2008) 

15. Takamizawa, K., Nishizeki, T., Saito, N.: Linear-time computability of combinatorial 
problems on series-parallel graphs. Journal of the ACM 29(3), 623-641 (1982) 


21 



